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C/3 \ Abstract 

We consider the Hirota equation (the discrete analog of the generalized Toda 
system) over a finite field. We present the general algebro-geometric method of 
construction of solutions of the equation. As an example we construct analogs 
of the multisoliton solutions for which the wave functions and the r-function 
can be found using rational functions. Within the class of multisoliton solutions 
^-j. , we isolate generalized breather-type solutions which have no direct counterparts 

O ' in the complex field case. 
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1 Introduction 

>- ! Cellular automata are mathematical models of physical systems in which space and 

time variables are discrete, and physical quantities take only finite number of values 



|34fl . In spite of simple formulation they are capable to describe wide variety of phe- 
nomena, for example traffic flow, immune systems, flow through porous media, fluid 
dynamics, ferromagnetism. Due to their completely discrete nature cellular automata 
are naturally suitable for computer simulations. However in this field there are not so 
many exact analytical results providing solutions with a given global behavior. The 
goal of this paper is to present a general method of construction of solutions to the 
cellular automaton associated with the Hirota equation. 



The Hirota equation fll8fl , the integrable discretization of the generalized Toda 
system [^7]], is one of the most important soliton equations. Its various limits give 
rise to variety of integrable equations, moreover it is the basic system for studying 
solvable quantum models pi] ]. 
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Like many other integrable systems the Hirota equation has simple geometric 
interpretation. In the paper we will use the following (geometric) form of the Hirota 
equation 

r m {n 1 ,n 2 ) r m (ni + l,n 2 + 1) = 

r m (ni + l,n 2 ) r m {n u n 2 + 1) - r m _ x {n x + l,n a ) W^n^na + 1). (1) 

One can describe it as the equation governing the so called Laplace sequence of two 
dimensional lattices of planar quadrilaterals |0 [[(J. This interpretation can be embed- 
ded into a more general theory of multidimensional lattices of planar quadrilaterals 
and their transformations [|]|, [15], [14], |l2j. It was noticed in || [£| that geometric 



constructions in the integrable discrete geometry (in particular, those leading to the 
Hirota equation) should work also on the level of finite geometries |L9fl, i.e., geometries 
over finite fields |Q. This observation has been developed in the present paper. 

The question of construction of integrable systems with solutions taking values 
in a discrete set (soliton cellular automata) is not new and it was undertaken in a 
number of papers, see for example [5|, || |33|, In particular, in paper g other 



Hirota equation (equivalent to the discrete sine-Gordon equation [17|) is investigated 
in the context of finite fields. 

In the paper we present general method of finding solutions of the Hirota equa- 
tion (|ID over finite fields by using algebro-geometric methods, standard in a complex 
domain in the soliton theory J22], ||]. We change however the field of definition of the 
underlying algebraic curves from the complex numbers to a finite field (see also earlier 
algebro-geometric papers ||, |32| where such a possibility was considered). 



It turns out that algebraic geometry over finite fields has become recently very 
important in practical use, especially in modern approaches to public key cryptog- 
raphy |20|] and in the theory of error correcting codes ||3l]]. With respect to the last 
application we would like to mention paper f[30|| where dynamics of the finite Toda 
molecule (a reduction of the Hirota equation) over finite fields was studied from the 
point of view of a decoding algorithm. 

We do not present here direct connection of the objects of the paper with the 
integrable discrete geometry over finite fields. This connection becomes clear when 
approach to the Hirota equation presented here (see also [22|]) is compared with results 
of U|, [11], ft^l , where methods of algebraic geometry over the field of complex numbers 
have been applied to construct integrable geometric lattices. 

The layout of the paper is as follows. In Section [2| we present the general algebro- 
geometric scheme for construction of solutions of the Hirota equation. Section |3| is 
devoted to construction of multisoliton solution on an algebraic curve starting from 
the vacuum solution. Finally, in Section |] we give in explicit form solutions of the 
Hirota equation for the background algebraic curve being the projective line. In 
particular, we present the mechanism (based on the action of the Galois group) of 
generation of generalized breather-type solutions and we discuss periodicity of the 
solutions. 
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2 Solutions of the Hirota equation from algebraic 
curves over finite fields 



This Section is motivated by algebro-geometric (over the complex field) approach to 
the Hirota equation (in a different form) |^2| and by papers [l], [TTJ, Q on application 
of algebro-geometric methods to integrable discrete geometry. It turns out that basic 
ideas of the algebro-geometric approach to soliton theory can be transferred to the 
level of integrable systems in finite fields. The notions and results form the theory of 



algebraic curves over finite fields which we use here can be found in |3T] . 

Consider an algebraic projective curve, absolutely irreducible, nonsingular, of 
genus g, defined over the finite field K = ¥ q with q elements, where q is a power 
of a prime integer p. By Cr we denote the set of K-rational points of the curve. By 
K denote the algebraic closure of K, i.e., K = IJfca^V? an d by denote the cor- 
responding (infinite) set of K-rational points of the curve. The action of the Galois 
group G(K/K) (of automorphisms of K which are identity on K) extends naturally 
to action on C^. 

Let us choose: 

1. two pairs of points a*, h G Ck, i = 1,2, 

2. N points c a G Cg-, a = 1, . . . ,N, which satisfy the following K-rationality condition 

Vct G G(K/K), a(c a ) = c a ,, 

3. N pairs of points dp,ep G Cg, (3 = 1,...,N, which satisfy the following Ir- 
rationality condition 

Va G G(K/K) : a({d p , e^}) = {dp, ep}, (2) 

4. g points / 7 G C^, 7 = 1, . . . , g, which satisfy the following K-rationality condition 

Va G G(K/K), a(/ 7 ) = f Y , 

5. the infinity point G Ck- 

Remark. We consider here only the generic case and assume that all the points used in 
the construction are generic and distinct. In particular, genericity assumption implies 
that the divisor D = Yl~ /= i f-y i s non-special. 

Remark. It is enough the check the K-rationality conditions in any extension field 
L D K of rationality of all the points used in the construction. 

Definition 1. Fix K-rational local parameters t{ at 6j, i — 1,2. For any integers 
nx,n2,m G Z define the wave function i/ji }7n (ni, ng) as a rational function with the 
following properties 

1. it has pole of the order at most n\ + m + 1 at 61 and the pole of order at most 
112 — m at 62, 
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2. its first nontrivial coefficient of the expansion in t\ at b\ is normalized to one, 

3. it has zeros of order at least n\ at ai, and of order at least n 2 at a 2 , 

4. it has at most simple poles at points c a , a = 1, . . . , N, 

5. it has zero at least of the first order at the infinity point hoo, 

6. it has at most simple poles at points / 7 , 7 = 1, . . . , g, 

7. it satisfies N constraints 

tpi,m{ni,n 2 ){d l3 ) = Tpi, m {ni,n 2 ){ep), p = l,...,N. (3) 

For the same set of points we define the wave function ip 2 ,m{ n i,n 2 ) as a rational 
function which differs from the function ipi, m (ni,Ti 2 ) only in the properties 1 and 2 
I2. it has pole of the order at most n\ + m at bi and it has pole of the order at most 
n 2 — m + 1 at b 2 , 

22. its first nontrivial coefficient of the expansion in t 2 at b 2 is normalized to one. 

Remark. The functions ipi,m(ni, n 2) are K-rational, which follows from K-rationality 
conditions of sets of points in their definition. 

Remark. As usual, zero (pole) of a negative order means pole (zero) of the corre- 
sponding positive order. Correspondingly one should exchange the expressions "at 
most" and "at least" in front of the orders of poles and zeros. 

Proposition 1. The wave functions ^, m (ni, n 2 ) are unique. 

Proof. We show this for the first function. By the Riemann-Roch theorem the di- 
mension (over K) of the divisor 

N g 

+ '^2 fj - hoo - rixai - n 2 a 2 + (m + 1 + m)6i + (n 2 - m)b 2 

a=l 7=1 

is equal to N + 1. Under the genericity assumption the N constraints (H) and the 
normalization condition at b\ remove the freedom. □ 

Corollary 2. In the next section we show, starting from the wave functions for N = 
one can construct the functions for arbitrary N. 

Remark. To make subsequent formulas more transparent from now on we will skip 
frequently the dependence on the parameters ri\,n 2 . 

In the generic case, which we assume in the sequel, when the order of the pole 
of ipi, m at 61 is rii + m + 1 and the order of the pole of ip 2 , m at b 2 is n\ — m + 1 
denote by Qi 2 ,m(n>i, n 2 ) the first nontrivial coefficient of expansion of ifti jm at b 2 , and 
by Q2i,m( n i> n 2) the first nontrivial coefficient of expansion of ijj 2)m at 61, i.e., 

1 1 

■01,m = n 2 - m (Ql2,m + ■■■), ^2,m = ni+m ($21,™ + •••)■ 

l 2 l 1 
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Remark. The functions Q\2, m an d Q2i,m take values in the field IK of the definition of 
the curve. 

Denote by Tj the operator of translation in the variable n,, i = 1,2, for example 
Proposition 3. The function ipi tm satisfies equations 

T 2 1pl, m - 1pl,m = (T 2 Ql2,m)lp2,m, (4) 

i>l,m+l - Tli>l,m = — ^l,m, (5) 

V12,m 

fpl,m-l = -^—^2,m- (6) 
V21,m 

T/ze analogous system for ip2,m is obtained by exchanging indices 1 with 2 and reversing 
the shift in the discrete variable m 

Tl1p2,m ~ 1p2,m = (TlQ21,m)lpl,m, (7) 

^2,m-l - T 2 ^2,m = — |j — 1p2,m, (8) 

V21,m 

^2,m+l = ^ V'l.m- (9) 
Wl2,m 

Proof. To prove the first equation (|]) notice that the left hand side has all properties 
of the function ip 2 ,m except of the normalization and must be therefore proportional 
to ^2,771- The coefficient of proportionality can be fixed comparing expansions at &2- 
Other equation can be proved in the same way. □ 

Fix K-rational local parameters £j at Oj, i — 1, 2. In the generic case when the 

order of ^ m at a; is rii by Pi, m (^i, ^2) denote first non-trivial coefficients of the 
expansion of ^i,m at a;, i.e., 

1pi,m ti (pi t m ■ ■ • ) • 

Similarly, in the generic case when the order of ipi jm at aj, j 7^ 2, is rij by Xij,m( n ii ^2) 
denote the first non-trivial coefficients of the expansion of ipi tm at a.j, i.e., 

1pi,m tj (Xij,m ~l~ • • • ) j i 7^ 3- 

Proposition 4. There exists a "K-valued potential (the r -function) defined (up to a 
constant) by formulas 

T\T m = p 1<m r m , (10) 

T2T m = P2,m,T m , (11) 
r m+l — ( — l) ni+n2 Ql2,mT"m- (12) 



5 



Proof. The first terms in expansion of equations ([|) and (|7|) at a\ give 

2~2Pl,m — Pl,m — {T2Q12 

— X21,m — {T\Q 21, m) Pi 1 

which combined together give 

T 2 Pl, m = (1 - (T 2 Qi2 im )(Ti<2 2 l,m))Pl,m- (13) 

Similarly, but changing the expansion point to a 2 we obtain 

T\p2,rn — (1 — (T 2 Q 12 t m)(TiQ 21 , m )) P2,m- 

Both equations imply 

{T2Pl,m)p2,m = (^lP2,m)pi,m, 

which is the compatibility condition of equations (|T0|) and ([11]). 
Expansion of equation (|5|) at ai gives 

Pl,m+1 = 7^ Pl,m, (14) 

V12,m 

which is the compatibility condition of equations ( |10| ) and flT2|) . Finally, by comparing 
equations (||) and @ we obtain 

,m+l f^. 1 

(15) 

V12,m 

which, combined with the following consequence of expansion of (||) at a 2 

T2Q21,m 

P2,m-1 — J\ P2,m, 

V21,m 

gives the compatibility condition of equations (|ll]) and (|12D . □ 



Corollary 5. Equation (|i3|) written in terms of the T-function reads 

T m TiT 2 T m = TiT m T 2 r m - Tir m _i T 2 r m+ i, (16) 
which is the Hirota equation ]7^/ . 



Corollary 6. Compatibility condition of equations (|T3|) and (0) written in terms of 
the function Qi2, m reads 

TiT 2 Ql2,m _ TlQl2,m _ T{IjQl2^n _ 7 2 <5l2,m+l iX?) 
^2<5l2,m Ql2,m 7lQl2,m-l Ql2 
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3 Construction of solutions using vacuum functions 

Results of this Section were motivated by papers on the fundamental transformation 



of quadrilateral lattices in a vectorial formulation |^6|, 15, |25| and on the algebro 



geometric interpretation of this transformation [|11], [12| . 

In the case N = let us add superscript to all functions defined above and call 
them the vacuum functions. The functions for arbitrary N can be constructed with 
the help of N new functions, which we define below. 

Definition 2. Fix local parameters t a at c a , a = 1, . . . , N. For any a define the 
function 0° m by the following conditions: 

1. it has pole of the order at most »i+m at bi and the pole of order at most n 2 — m 
at 6 2 , 

2 it has zeros of order at least ri\ at a\, and of order at least n 2 at a 2 , 

3. it has at most simple pole at the point c a and the first nontrivial coefficient of the 
expansion in t a at c a is normalized to one, 

4. it has zero at least of the first order at the infinity point hoc, 

5. it has at most simple poles at points / 7 , 7 = 1, . . . , g. 

Remark. The function 4>° a m is unique but usually it is not K-rational. 

Lemma 7. Denote by ifj^ m (d, e), i = 1, 2, the column with N entries of the form 

[< m (d, e)] p = < m (^) - < m ( e/3 ), /3 = 1,...,N, 
denote by <j>° A m the row with N entries 

[Am] = C» a=l,...,N, 
and denote by (f> Am (d, e) the N x N matrix whose element in row a and column (3 is 

[KM e)] ap = 4>l m (dp) - < m ( e/9 ), a,P = l,...,N. 
Then the wave functions i[)i }Tn read 

Am = V?m - <P%J<P% m (d, e)]-VSJd, e). (18) 

Proof. Denote the right hand side of equation (|i~8D by ipi m . By building the column 

tpi m (d,e) with iV entries of the form i$ m (dp) — tpi m {&p) we can easily show that 

ipf m (d,e) = 0. This demonstrates that the function ip^ m satisfies constraints (|||). 

One can check check that ip® m satisfies also other properties which define uniquely 
the function ipi m . □ 
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In the generic case denote by Hf a m the first nontrivial coefficient of expansion of 
the function in the uniformization parameter U at fcj, for example, 



i>° = 

a,m ,n 



1 



i+m \^l,a,m + • • • ) • 



Lemma 8. The corresponding expressions for Qij, m , i 7^ j> and for pi :m read 

Qij,m = ^-^J^elWe), i±h (19) 
p hm = p° m (l + (T^ A J[0 o Am (d,e)]-V° m (d,e)), (20) 



where H^ A m is the row with N entries Hf 



Proof. Because equation fll9|) can be obtained by expansion of formula fll8|) at bj, 
only equation ( f20| ) needs an explanation. Notice first equation 

Ti4Ln ~ Cm = (TiH°Jl£ m , (21) 







a, m 



which can be shown in the same way like equations of Proposition |3|. Denote by F- 
the first nontrivial coefficient of expansion of 0° at a*, for example 

fh° - — (F° 4- 1 

Ta,m 7ni \ l,o,m '"'')' 

Expansion of equation fll!f ) at a$ gives 

ri,m\ i i,a,m) i,a,mi 

which concludes the proof. □ 

We will use the following result, which can be proven by induction with respect 
to the dimension of the vector space V. 

Lemma 9. Given u e V and v* G V* ; if ly is the identity operator on V then 

det(l v + u <g> «*) = 1 + (v*, it). 

Proposition 10. Using the above notation the T-function can be constructed by the 
following formula 

T m = T m detA m (d,e). (22) 



Proof. Notice that equation fl21|) implies 

[<t>\ m {d, e)r x T4\ m {d, e) = V + ([^ m (d,e)]-V; m (d,e)) ® (T^J, 
which gives, by Lemma |9|, 

^Ittjd,e) = 1 + & H ^rn)l<f>\M eT^Jd, e). (23) 



Comparing equation fl23| ) with equation (|20| ) and taking into account equations (p~0|)- 
(|TlD we obtain 

det r^ m (d, e) _ X t r m /T t T° 

det^ m (d,e) r m /r° ' 1 j 

Notice the following equation 

H° 

Vl2,m 

which can be shown in the same way like equations of Proposition ||. It implies that 
[0A, m (d,e)]-V°A, m+ i(d,e) = Igs - ([^ m (d,c)]-V? iro (d,c)) ® (# 2 ° A J, 

Vl2,m 



which gives 

det0^ m+1 (d, e) i M-i ,o 



Hpt r r. ' ' = 1 - y^H^JcPXJd, e)}-^lM e). (25) 



Comparing equation (|25|) with equation fll9|) and taking into account equation (|12| ) 
we obtain 

det <&, m+ i(d, e) = r m+1 /r° +1 

det^ m (d,e) r m /r° ' 1 Dj 

which, together with equation (^4j), concludes the proof. □ 



Corollary 11. Notice that equation fl22 ) is wa/id up £o a (nonessential) change of the 
initial value of the r -function, which is due to introduction of the integration constant 
from formulas (|24]) and fl26| ) . 

Corollary 12. Starting with K-valued function r° t and £/ze local parameters t a at c a 
chosen in a consistent way with the action of the Galois group G(K/K) on we 
obtain K-valued function r m . 



4 Multisoliton solutions 

We present here explicit formulas for the vacuum functions in the simplest case g = 
0. Then with the help of these expressions we present some examples of iV-soliton 
solutions. In constructing the vacuum functions we will use the standard parameter 
t on the projective line P(K) and we put = oo. 
Explicit forms of the wave functions read 



n, 



o 1 (t-aiY^t-atfvQn 



(t- b 1 ) n ^+ 1+m (6i - ai) ni (6i - a 2 )" 2 (t - 6 2 ) n2_m ' 



"2 



o 1 (&2-&i) ni+m (*-ai) wl (*-«2; 

^ 2 ' m (t - 6 2 ) na+1 - m (t - 6i)" 1+m (& 2 - ai) ni (6a - a 2 ) n2 ' 
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which gives formulas for the functions Q\ 2 m and Q 2 



21,771 



Ql2,m 



and for the functions pi >m and P2, 



Pl,m 





/ 1 \ T7 n — IT) 
1 y^A '"• 




(62 - 


ai) H & 2 - 


- a2)™ 2 


(&2 






/ 7 

(6l " 


\ IF, / T 

ai)"i(6i - 


- a 2 ) n2 




I 1 IL\-\-tlL 




/ 7 


ai) ni (6i - 


- a2) n2 




— 5 2 ) n2_ni + 1 


—2m 


(62 " 


ai) ni (& 2 - 


- a 2 ) n2 


and 


P2,m 












(_l)m 


(ai 


- 02 




jn 2 — 771 


(ai - 




(61 


- 


712 (at - b 2 


J 712— 771 ' 




(_l)n 2 


(62 


-61) 


ni+m (a 2 - 


ai) ni 


(«2 " 




(^2 


-61: 


"i+ m (6 2 - 


ai) ni ' 



Explicit form of the vacuum r-function reads 

( ^_l^lni(m-l)+n 2 (n2-i)+m(m+l)]/2 ^ _ q^™!"!!^ _ ^ (n 2 -m)(m+rn) 



(ai - 6i) ni ( ni+m )(a 2 - 6 2 ) n 2( n 2-m) _ a 2 )n 2 (ni+m)( ai _ /j 2 )m(n 2 -m) ' 



The functions 0" ^, a = 1, . . . , iV have the form 

= 1 (t ~ ai) ni (t - a 2 ) n2 (c a - 61)™+^ - 6 2 )" 2 - m 
,m * - c a (c a - ax)"i (ca - a 2 )" 2 (t - - b 2 ) n *~ m ' 

and can be used, due to Proposition to construct the r-function for arbitrary N. 

Let L = ¥ q e cKbea field of rationality of all the points used in the construction. 
Recall [23|] that if q = p k then the Galois group G(L/K) is the cyclic group of order 
£ and is generated by ap, where ap is the Frobenius automorphism of L defined as 
(Tp(a) = aP . Therefore possible systems of points c a and pairs {dp, e@} can be grouped 
into K-rational clusters (orbits of the group G(L/K)) of the lengths being divisors 
of I. In the standard nomenclature the clusters of length one correspond to solitons, 
and clusters of length two give rise to breathers (the position of poles c a of the wave 
functions must be symmetric with respect to the complex conjugation). In finite fields 
we encounter new types of solutions (without direct analogs in the complex field case) 
which come from clusters of lengths greater then two. The analog of the breather 
solution will be presented in Example |l|. Let us call the iV-soliton solution of order 
I the K-rational iV-soliton solution with parameters in extension of IK of order £. In 
this terminology the standard iV-soliton solutions are of order one, while the breather 
solution is a 2-soliton solution of order two. A 3-soliton solution of order three and a 
2-soliton solution of order four are presented in Examples || and ||. We remark that 
the above terminology is not completely distinctive. 

Notice that the variables n\, n 2 and m enter exponentially in the functions r^ and 
4>° a m This implies that the r-function is periodic in rii, n 2 and m with the periods 
being divisors of q l — 1, which is the order of the cyclic multiplicative group L*. 
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Finally, we present the examples. For any example we describe the fields K and 
L giving first the numbers q = p k and £ and then writing down the polynomial w(x) 
over F p used to construct multiplication in the field L. We represent elements of L 
as elements of the vector space F* 1 . Then we give the points used in the construction 
of the solution of the Hirota equation presenting also the action of the Galois group 
G(L/K) on them. 

Example 1. A breather solution of the Hirota equation in F 5 . Parameters of the 
solution take values in extension F 5 2 of F 5 by the polynomial w(x) = x 2 + x + 1. 
The corresponding Galois group reads G(F 5 2/F 5 ) = {id, a}, where a 2 = id. The 
parameters of the solution are chosen as follows 

ai = (00), a 2 = (02), h = (01), b 2 = (04) 

ci = (10), c 2 = a( Cl ) = (44), 

da = (11), d 2 = aid!) = (40), 

ei = (13),e 2 = a( ei ) = (42). 

This solution is presented in Figures [I] and ||. The elements of F 5 are represented 
by: ■ - (00), ■ - (01), (02), ■ - (03), ■ - (04). 

The periods in variables m, n 2 and m are 12, 24 and 24, correspondingly. Notice 
that figure for m = 5 can be obtained from the figure for m = 1 by a shift in n 2 by 
four. 

Example 2. A 3-soliton solution of order 3 of the Hirota equation in F 5 . This 
solution is the first one in the generalized breather class, a new type not present in 
the standard complex field case. Parameters of the solution take values in extension 
F 5 3 of F 5 by the polynomial w(x) = x 3 + x 2 + 1. The corresponding Galois group 
reads G(F 5 3/F 5 ) = {id, a, a 2 }, where a 3 = id. The parameters of the solution are 
chosen as follows: 

ai = (004), a 2 = (003), h = (002), b 2 = (001) 

ci = (022), c 2 = (7(d) = (120), c 3 = a 2 ( Cl ) = (412), 

di = (020), d 2 = a(di) = (123), d 3 = a 2 ^) = (410), 

e x = (010), e 2 = a{e x ) = (314), e 3 = a 2 ( ei ) = (230). 

This solution is presented in Figures |] and |]. The elements of F 5 are represented 
like in Example [I]. From the figures one can deduce that the periods in all variables 
ni, n 2 and m are maximal and equal to 124. 

Example 3. The goal of this example is to present a solution in a "small" field F 4 
obtained from parameters taking values in a relatively bigger field F 256 . The field 
F 256 is chosen as extension of F 2 by the polynomial w(x) = x 8 + x 6 + x 5 + x + 1. 
The elements (00000000), (00000001), (11110000) and (11110001) of F 256 = F 2 s form 
a subfield isomorphic to F 4 . The Galois group is generated by a — a\, and reads 
C(F 2 5 6 /F 4 ) = {id, a, a 2 , a 3 }, where a 4 = id. The parameters of the solution are as 
follows: 

ai = (00000000), a 2 = (00000001), b x = (11110000), b 2 = (11110001) 
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Figure 1: A breather solution of the Hirota equation in F 5 described in Example [Q; 
rii range from to 28 (directed to the right), 712 range from to 28 (directed up), 
m = 1. 



Figure 2: A breather solution as in Figure [I] for m = 2 and m = 5. 
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Figure 3: A breather of order 3 solution of the Hirota equation in F 5 described in 
Example 0; n-y range from to 8 (directed to the right) , n 2 range from to 8 (directed 
up), m = 0, 4 and 62. 

ci = (00010010), c 2 = ct(ci) = (11100011), 

di = (00001010), d 2 = a{d x ) = (00001001), 

ei = a 2 ^) = (00011000), e 2 = <x 3 (di) = (11101010). 

Here the points C\ and c 2 are chosen from a subfield of F 256 isomorphic to F 16 and 
form a cluster of length two. The points ci, c 2 , di and <i 2 form a cluster of length 
four in a way compatible with the F 4 -rationality condition @. We obtain therefore a 
2-soliton solution of order four, which also has no direct counterpart in the complex 
field case. 

The solution is presented in Figure |5| for m = 0; The elements of F 4 are represented 
by: ■ - (00000000), ■ - (00000001), ■ - (11110000), - (11110001). The period 
in all variables is the same and equals 51. 

5 Conclusion and remarks 

In the paper, motivated by recent developments of integrable discrete geometry, we 
presented algebro-geometric method of construction of solutions of the Hirota equa- 
tion over finite fields. It turns out that the main ideas used for integrable systems 
over the field of complex numbers, e.g., application of the Riemann-Roch theorem, 
can be transferred to the level of finite fields without essential modifications. We 
would like to stress that although finite fields consist of finite number of elements the 
corresponding algebraic curves have infinite number of points (taking into account 
the algebraic completion of the field) which gives rise to infinite families of solutions 
of the equation. 

We have presented examples of pure (for the simplest algebraic curve being the 
projective line) multisoliton solutions of the Hirota equation. Less trivial examples 
which use techniques on Jacobians of algebraic curves will be the subject of another 
paper @. 
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Figure 4: A "global" view of the first picture (m = 0) of the generalized breather 
solution presented in Figure |3|; n\ and n 2 range from to 129. 
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Figure 5: A 2-soliton solution of order four of the Hirota equation in F 4 ; n\ and 712 
range from to 59, m = 0. 
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